Edge pinch instability of liquid metal sheet in a transverse high-frequency ac magnetic 

field 
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We analyze the linear stability of the edge of a thin liquid metal layer subject to a transverse 
high-frequency ac magnetic field. The layer is treated as a perfectly conducting liquid sheet that 
allows us to solve the problem analytically for both a semi-infinite geometry with a straight edge 
and a thin disk of finite radius. It is shown that the long-wave perturbations of a straight edge 
are monotonically unstable when the wave number exceeds the critical value kc — -Fb/(7^o), which 
is determined by the linear density of the electromagnetic force -Fo acting on the edge, the surface 
tension 7, and the effective arclength of edge thickness Iq. Perturbations with wavelength shorter 
than critical are stabilized by the surface tension, whereas the growth rate of long- wave perturbations 
reduces as ~ for — > 0. Thus, there is the fastest growing perturbation with the wave number 
fcmax ~ 2/3fcc. When the layer is arranged vertically, long- wave perturbations are stabilized by the 
gravity, and the critical perturbation is characterized by the capillary wave number kc — \/ gp/j, 
where g is the acceleration due to gravity and p is the density of metal. In this case, the critical linear 
density of electromagnetic force is -Fb,c = 2fccZo7, which corresponds to the critical current amplitude 
Io,c = 4-y/ TrkcloLj/fio when the magnetic field is generated by a straight wire at the distance L 
directly above the edge. By applying the general approach developed for the semi-infinite sheet, we 
find that a circular disk of radius Ro placed in a transverse uniform high-frequency ac magnetic field 
with the induction amplitude -Bo becomes linearly unstable with respect to exponentially growing 
perturbation with the azimuthal wave number m — 2 when the magnetic Bond number exceeds 
Bmc = -Bo-Ro/(2a'o^o7) ~ 37r. For Bm > Bmc, the wave number of the fastest growing perturbation 
is mmax ~ [2Bm/(37r)] . These theoretical results agree well with the experimental observations. 



PACS numbers: 47.20.Ma, 47.65.-d, 47.10.A- 



I. INTRODUCTION 

In several induction melting processes, such as the cold 
crucible or electromagnetic levitation, liquid metal with 
a free surface is subject to ac magnetic fields that may 
cause considerable deformations of Hquid metal resulting 
from the electromagnetic forces due to the eddy currents, 
which are often confined in a thin skin layer beneath the 
surface P|. It has been observed that the free surface 
sometimes may become strongly asymmetric and even 
irregular when a sufficiently strong magnetic field is ap- 
plied HESIil- 

Most of theoretical studies of the effect of ac magnetic 
field on the stability liquid metal surfaces have been con- 
cerned with fiat surfaces subject to tangential uniform 
magnetic field. McHale and Melcher were the first to 
show that the time-averaged electromagnetic force has a 
destabilizing effect giving rise to traveling waves on the 
surface of Hquid metal. Although the theoretical insta- 
bility threshold is in good agreement with experimental 
results, the predicted growth rates are too small com- 
pared to the experimental observations. Note that such 
small growth rates are typical for the electromagnetic in- 
stabilities caused by the currents induced by motion of 
conducting media in ac magnetic fields Q. This simple 
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model was revisited by a number of authors using vari- 
ous approximations. First, Garnier and Moreau found 
that ac magnetic field has only a stabiHzing effect on the 
surface waves when the currents induced by metal fiow 
are neglected. Deepak and Evans Q took into account 
the motion of a surface but not the associated fiow in the 
liquid, although both have a comparable effect, and they 
concluded that ac magnetic field can, however, give rise 
to unstable traveling surface waves. Stability of a fiat 
metal layer suspended by means of a uniform magnetic 
field, as in the experiment of Hull et al. J?!, has been 
studied by Ramos and Castellanos who analyzed the 
effect of the viscosity on Rayleigh- Taylor type instabil- 
ity, which is unavoidable in this system. Fautrelle and 
Sneyd ^ used a more elaborate model, taking into ac- 
count not only the time-averaged but also the oscillating 
part of the electromagnetic force, which results in much 
stronger parametric instabilities when the frequency of 
surface waves is sufficiently close to the multiple of the 
electromagnetic force frequency. Note that this simple 
model of a fiat surface with tangential magnetic field 
leads only to traveling but not stationary wave instabil- 
ities, which require consideration of nonplanar surfaces 
in nonuniform magnetic fields. A stability analysis was 
performed by Karcher and Mohring |i.l2] to describe the 
experimental observations of static surface instabilities 
by Mohring et al. However, drastic simplifications 

were made in the latter analysis. First, the authors used 
the mirror image method to find the magnetic field distri- 
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x(z,t)=^t(t)cos(kz) 



summarized in Sec. IV. 



II. PERFECTLY CONDUCTING LIQUID 
SHEET 



Mathematical model 




Figure 1: Sketch to the formulation of the problem with the 
magnetic flux lines (dashed) bending around perturbed edge 
x(z, t) — x{t) cos(fcz). 



bution at the end of annular gap filled with liquid metal, 
however this method is applicable only to simple geome- 
tries, such as half-space or a sphere [1^. Second, they 
neglect the effect of surface perturbation on the mag- 
netic field distribution, although the coupling between 
both constitutes the basic mechanism of this instability. 

In this work, we propose a simple theoretical model 
to describe such static surface instabilities. The model 
consists of a flat liquid metal layer in a transverse ac 
magnetic field. The ac frequency is assumed to be high 
so that the magnetic field is effectively expelled from the 
layer by the skin effect. The layer is assumed to be thin 
so that it can be regarded as a liquid perfectly conducting 
sheet. We start with the linear stability analysis of the 
straight edge of a semi- infinite liquid sheet, which allows 
us to work out a general approach, which is applied later 
to the liquid layer in the magnetic field of a straight wire 
parallel to the edge and to a thin circular liquid disk 
in a uniform transverse magnetic field. We describe a 
pinch-type instability of the edge of a liquid metal sheet 
with the following mechanism. The magnetic field bends 
around the edge of a perfectly conducting sheet, giving 
rise to the magnetic pressure on the edge, which tries to 
compress the sheet laterally. In the case of a straight 
edge, the uniform magnetic pressure along the edge is 
balanced by a constant pressure in the sheet. A wavelike 
perturbation of the edge causes the magnetic fiux fines to 
diverge at wave crests and convergence at troughs. This 
redistribution is because the magnetic fiux lines along the 
sheet are perpendicular to the electric current lines that 
are directed along the edge and, thus, the magnetic fiux 
lines have to be perpendicular to the latter. As the re- 
sult, the magnetic pressure is reduced at the crests and 
increased at the troughs, which enhances the perturba- 
tion. 

The paper is organized as follows. In Sec. II, we con- 
sider the general model of a semi-infinite perfectly con- 
ducting liquid sheet with a straight edge, and compare 
with the experiment. In Sec. Ill, the linear stability 
problem for a thin disk in a transverse magnetic field is 
solved and compared with experiment. The results are 



Consider a thin layer of liquid metal submitted to a 
transverse ac magnetic field B. We assume that the 
layer is semi-infinite and lies on the right-hand side of 
x-z plane, so that the unperturbed edge of the layer co- 
incides with z axis, as shown in Fig. ^ Field frequency 
is assumed to be high, so that the layer is effectively 
impermeable to the magnetic field because of the skin 
effect. In addition, the layer is assumed to be a thin 
sheet in a static equilibrium state supported by a fiat 
horizontal non-wetting surface or constrained in the gap 
between two parallel walls. Note that this model repre- 
sents a special case of a thin superconductor film (see, 
e.g., ^3)- The magnetic field B in the free space around 
the sheet can be described by the scalar magnetic poten- 
tial ^. Then the solenoidity constraint V • B = results 



V^* = 0. 



(1) 



The impermeability condition at both sides of the sheet 
is 



(n-B 



\y=±Q;x>0 — 



0, 



(2) 



y=±Q; x>0 



where n is the surface normal vector. First, we focus 
on the distribution of the magnetic field in the vicinity 
of the edge, which can conveniently be described in the 
cyHndrical coordinates with the z axis coinciding with the 
edge and the polar angle (p measured from the 
illustrated in Fig. The solution for the unperturbed 
potential in the vicinity of straight edge satisfying condi- 
tion (El is El 



•9o{r,(p) = CoVrcos{tp/2), 



(3) 



where Cq is an unknown constant. According to sim- 
ple dimensional arguments, determining Cq requires an 
external length scale that, however, is missing in this 
simple model. Thus, Co can be determined for a strip of 
finite width, a magnetic field generated by a straight wire 
placed at some distance parallel to the edge or a circular 
disk of finite size, that will be done in the following sec- 
tion. But first, we develop a general approach without 
specifying Cq. 

Suppose that there is a perturbation of the edge posi- 
tion X = xi{z,t) — x{t) cos{kz) with a small, generally 
time-dependent amplitude x{t) and the wave number k 
along the z axis. This perturbation gives rise to the po- 
tential perturbation that can be presented as 

^'(r, ip, z) = ^'o(r, p) + £*i(r, z) + • • • , 
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where ^'i is a perturbation with small amplitude e. To 
relate the perturbation of a potential to that of the edge, 
we need an additional condition at the edge, which is 
derived as follows. For the surface current with density 
J, we have /ioJ — n x B where fj,o is the perme- 
ability of vacuum. According to this relation, the mag- 
netic field along the sheet is perpendicular to the cur- 
rent. Consequently, the magnetic field has be to per- 
pendicular to the edge because the current is fiowing 
along the latter. Thus along the edge L, we obtain 
{t ■ B)\j^ = 95'/9r|^ = 0, where r is the unit vector 
tangential to the edge. This condition in turn implies 
that^"!^ — const, where const = may be chosen be- 
cause the potential is defined up to an additive constant. 
Applying this condition at the perturbed edge, we obtain 
up to the first-order terms in the perturbation amplitude 



which results in 



9*0 

dx 



■Xi 



0, 



9*0 
dx 



■Xl 



r^O; ip=0 



r-^0; 1^=0 

_ Co Xl 
2 



(4) 

Note that the base potential cannot formally be expanded 
in power series directly at the edge because the necessary 
derivative is singular there. To circumvent this, we take 
the expansion not exactly at the edge but define it as a 
limit when the expansion point approaches the edge. 

The potential perturbation that satisfies Eq. ^ is of 
the same form as the base field *i(r, Lp) — fi{r) cos{(p/2). 
When substituted into Eq. 1^, this leads to 

r dr \ dr I 4 



The solution 

-\CaXte-''-' 
base field is 



satisfying Eq. Q is fi{r) — 
and the full potential including the 



1 Xl 



-kr , 



cos(fcz) I COs((y9/2). 



(5) 

Note that the potential above, which is defined relative 
to the unperturbed edge, contains a singularity at the 
unperturbed edge. This singularity can be removed by 
proceeding to the coordinates defined relative to the per- 
turbed edge as r = r' + exXi, and expanding the solution 
in terms oixi. Thus, we obtain up to the first-order terms 
in the perturbation amplitude 



^(r'-l-e^xi) 



dx 



Xl + e*! 



(6) 



bvr' 



1 



2r' 



1 



r=r 

kr 



cos{kz] 



where r' is the cylindrical radius relative to the perturbed 
edge. Having no singularity anymore, this solution sim- 
plifies in the vicinity of the edge to 

*L, , 




y 



y 



CoVr' cos((^/2). 



Figure 2: Scalar magnetic potential in the vicinity of a 
straight edge (a) and its perturbation amplitude with l/fc 
used as the length scale (b). 



where Cq = Cq [l -I- ^xk cos(fcz)] . Thus perturbation of 
the edge results just in a redefinition of the constant Co, 
which is now replaced by Co, whereas the distribution of 
the potential remains the same as for the straight edge 
obtained above. This is because a smoothly perturbed 
edge looks straight again when examined on a sufficiently 
small scale. The scalar magnetic potential in the vicinity 
of a straight edge and its perturbation amplitude defined 
by Eq. ^ are plotted in Fig. |2lwith l/fc used as the 
length scale. The corresponding magnetic fiux and cur- 
rent lines along the layer in the vicinity of the perturbed 
edge are shown in Fig. El As seen, the magnetic fiux 
lines diverge at wave crests and converge at troughs in 
order to remain perpendicular to the edge, as discussed 
above. This redistribution of the magnetic fiux lines at 
the perturbed edge is the principal physical mechanism 
behind the instability considered in this study. 

In the perfect conductor approximation, the electro- 
magnetic force due to an ac magnetic field reduces to an 
effective magnetic pressure acting on the surface of the 
layer with the time-averaged value Pm = |-Bo| /(4/io), 



4 




Figure 3: Magnetic flux and current lines along the layer in 
the vicinity of the perturbed edge. 



where Bq is the amplitude of an ac magnetic field. Note 
that part of the magnetic pressure, which oscillates with 
double ac frequency, is neglected here by assuming the 
frequency to be so high that inertia precludes any consid- 
erable reaction of the liquid. According to Eq. 10, the 
magnetic pressure increases toward the edge as ~ 1/r 
and, thus, it becomes singular at r = 0. Nevertheless, 
the integral force on the edge has a finite value. This 
is because the magnetic pressure at the edge of a layer 
of small but finite thickness do increases as ~ l/^Oj 
which, integrated over the thickness, results in a finite 
value independent of do. The force on the edge can be 
evaluated by integrating the Maxwell stress tensor over a 
small cylindrical surface S enclosing the edge, as shown 
in Fig. U that results in the first-order terms in 



lim 





Fi) 



2 



(B n)B 



rdip 



where 



-Fo = 



8^0 



(7) 



(8) 



and Fi = Fgkx cos{kz) are the base force and its pertur- 
bation, respectively. 

Further, we assume the sheet to be an inviscid liquid, 
and consider a small-amplitude potential fiow associated 
with the edge perturbation. Justification of this assump- 
tion, particularly at the threshold of monotonous insta- 
bility, is discussed at the end of the section. Then the 
linearized Euler equation applied to a potential velocity 
field V = V$, 



dv 



Vp^V [p 



~dt 



-P 



0, 



leads to the pressure distribution in the sheet p = Po — 
pd^/dt = Po + pi, where po is a constant base pressure 




Figure 4: Evaluation of the electromagnetic force on the 
edge by integration of the Maxwell stress tensor over a small 
cylindrical surface S enclosing the edge. 



and pi ~ —pd^/dt is the perturbation of pressure. Ve- 
locity potential <& is governed by the incompressibility 
constraint 'V ■ v — resulting in 



(9) 



We integrate the normal stress balance over the edge by 
assuming both the pressure and curvature to be constant 
because of the small thickness of the layer, which yields 



P\x=0 = 



1 

R 



F 
To' 



(10) 



where lo is the effective arclength of the edge thickness; 
7 is the surface tension and l/R denotes the curvature 
of the edge. For an unperturbed edge, we have po = 
7/i?o + Fo/lo- Then for the perturbation, the balance 
condition takes the form 



'dt 



x=0 



1_ 

Ri 



F\ 



(11) 



where 1/Ri ~ V^xi is the curvature perturbation of the 
edge. In addition, we have a kinematic constraint 



la:=0 



dx 



dxi 
'dt' 



(12) 



Now we search for the amplitude of edge perturbation of 
the form x{t) = xoe"^*, where A is, in general, a complex 
growth rate, whose real part has to be negative for the 
perturbation to be stable. The hydrodynamic potential 
is of the form 

$(x, z, t) = l>(x) cos(fc2;)e^*, 

which when substituted into Eq. Q leads to d^^/dz^ — 
= 0, whose solution decaying away from the edge is 
^{x) = <i>oe~'^^. The amplitude of the hydrodynamic po- 
tential is related to that of the edge perturbation by the 
kinematic constraint II12II $o — ^xoX/k. Finally, the nor- 
mal stress balance Ijllll yields the growth rate depending 
on the wave number 



(13) 



5 




which implies that long-wave perturbations with the 
wave numbers < k < kc = -Fo/(7^o) have positive 
growth rates and, thus, they are unstable, as illustrated 
in Fig. The stronger the electromagnetic force Fq 
on the edge, the shorter the critical wavelength. The 
waves that are shorter than the critical one are stabi- 
lized by the surface tension. Although long waves are 
always unstable, their growth rate reduces as ~ fc for 
A: ^ 0. Thus there is a perturbation with /cmax = f 
for wh ich the gro wth rate attains the maximum Amax — 

kmaxy/Fo/{3plo) (sCC Fig. EJ. 

Concerning the effect of neglected viscosity, simple 
physical arguments suggest that this instability can only 
be slowed down but not prevented by the viscosity. Note 
that the viscosity is inherently related to the fluid flow. 
But at the threshold of monotonous instability, where 
X{kc) = 0, the characteristic time of monotonous insta- 
bility tends to infinity. Consequently, there is no charac- 
teristic time and, thus, no characteristic velocity scale for 
the monotonous marginally stable mode which, therefore, 
cannot be affected by the viscosity. 




B. Comparison with experiment 



To compare our theor y w ith the experiment of Mohring 
et al. [JJj, similarly to we unfold the annular layer 
of InGaSn (Galinstan) melt used in experiment by ap- 
proximating it by a semi-infinite fiat perfectly conducting 
liquid sheet. Magnetic field is approximated by that of a 
straight wire lying at the distance L from the edge in the 
plane of the sheet. The distance L provides us with the 
length scale necessary to specify the constant Cq used in 
our model above. This constant follows from the com- 
plex potential of the magnetic field, which is obtained by 



Figure 6: Magnetic flux lines and isolines of the scalar mag- 
netic potential represented by the real and imaginary parts of 
the complex potential for a straight wire placed at a; = L = 1 
and a semi-inflnite sheet at a; < (a) and for the sheet of fl- 
nite height —H < x <Q transverse to a perfectly conducting 
wall At X = -H ^ -2 (b). 



the conformal mapping 



w{0 



■log 



2tt 



where C = x + iy, as Co = ^iqIo/ {tt^/L). The magnetic 
fiux lines and isolines of the scalar magnetic potential 
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Figure 8: Sketch of a thin perfectly conducting disk in axial 
ac magnetic field. 



Figure 7: Critical current amplitude vs the layer thickness. 
The circles correspond to the measured total current ampli- 



tude for the g ap l^c 



17 mm between the metal surface and 



the lower side of the inductor 11]. The theoretical curve cor- 
responds to the best fit with the effective gap with L ^ 45 mm 
for a semi-infinite sheet and L ~ 28 mm for a sheet of finite 
height H = 24 mm with a perfectly conducting bottom. 



represented by the real and imaginary part of this com- 
plex potential are shown in Fig. Eta) with L used as the 
length scale. In addition, we consider the gravity with 
the acceleration g = e^g directed downwards along the 
sheet normally to its edge. Then Eqs. ltT3|l and l|Hl with 
Co defined above result in 



X{k) = kJ- 



1 f f^olo 



p \8t:LqIo 



9 
k' 



(14) 



As is easy to see, the gravity stabilizes long-wave dis- 
turbances, whereas short waves are stabilized by the sur- 
face tension. Thus, the first unstable mode, defined 
by A(fcc) — 0, appears at the capillary wave number 

kc = ~ 0.296 mm^^ that corresponds to the crit- 

ical wavelength Ac = 27r/fcc « 21.2 mm, where p = 
6440 kq/m-^ and 7 = 0.718 N/m are the density and sur- 
face tension of Galinstan Note that in the exper- 
iment, the surface of liquid metal is covered by a layer 
of NaOH solution. Thus, perturbation of the hydrostatic 
pressure at the interface is determined by the density 
difference of GalnSn and NaOH. Assuming the latter to 
have the density of water, we find the critical wavelength 
Ac ~ 23 mm that coincides very well with that of the 
static surface deformation observed in the experiment. 
The critical electromagnetic force follows from Eq. Ijl4|l 
as -fo,c = \/27^o/fcc, which corresponds to the critical 
current amplitude 



Po 



(15) 



approximated by a half-circle, i.e., Iq — ■ndQ/2. In or- 
der to compare this result with the measured critical 
currents^3|) note that the coil used in the experiment 
consists of two horizontal layers, each of which contains 
five windings. Thus, the measured current has to be mul- 
tipHed by 10 to obtain the total current amplitude. Un- 
fortunately, the authors do not specify the coil dimen- 
sions but give only the gap width Lcxp = 17 mm between 
the metal surface and the lower side of the inductor which 
is not sufficient for comparison with our theory. There- 
fore we treat the distance L as a free parameter to fit 
the experimental results that yield L « 45 mm (see Fig. 

. Note that the model of a semi-infinite layer may not 
be very adequate for the given experiment with the layer 
extention H — 24 mm which is comparable to the gap 
width L, especially when the layer resides on a well con- 
ducting metal plate. Finite extention of the layer and 
the conducting bottom can partly be accounted for by a 
more sophisticated complex potential. 



w{0 



where the edge arclength Iq over the layer thickness do is 



which is plotted in Fig. Efb). This yields Co = 
Polo/ jTTv/i [1 + L/{2H)]j , resulting in L « 28 mm 
which is considerably closer to the corresponding exper- 
imental value. 

III. A THIN LIQUID DISK 

A. Analytical solution 

Now we will apply the approach developed in the pre- 
vious section to a thin circular liquid disk of radius Ro 
and fixed thickness do, which is subjected to a uniform 
axial ac magnetic field with an induction amplitude Bq, 
as shown in Fig. |H| The thickness c?o is assumed to be 
small relative to the radius of disk do ^ Ro, so that the 
disk may be regarded as a thin sheet. The magnetic field 
is sought in terms of the scalar magnetic potential ^' gov- 
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erned by Eq. (QJ, whereas the impermeabihty condition 
at the disk surface S takes the form 



dn 



0. 



(16) 



At large distances from the disk, the field is uniform and 
axial, which implies 



|r|- 



(r • Bo) = zBo, 



(17) 



where z is the axial distance from the disk. Solutions for 
both a circular and a slightly perturbed disk can be ob- 
tained analytically in the oblate spheroidal coordinates, 
which are related to the cylindrical ones by 



z = Rov£„ 

where < 77 < 1 and < <^ < 00 are the angular and 
radial spheroidal coordinates, respectively, as defined in 
[r^ . Equation 1^ for the unperturbed potential ^'o 
around a circular disk takes the form 



d_ 



(1 



9£o 
drj 



d_ 



1+e 



The impermeability condition Ijl6|l is 



9C 



= 0. 



0. (18) 



(19) 



The second boundary condition Ijl7|l suggests a solution 
of the form *o(?7, = vMO- This results in [l^ 



*o(?7,0 =^^o7?[l + earctan(e)]. 



where Co ~ 2RoBo/n, which is plotted in Fig. El with 
the corresponding magnetic flux lines. Note that in the 
vicinity of the edge this solution reduces to ^^o = CqT] + 
0(^^), which is equivalent to Eq. JSJ. 

Further, let us consider a perturbation of disk radius 
along the azimuthal angle 4> of the form 

Ri<j),t)^Ro + RTi<l),t), 

where RY^{(f>,t) = i?™(t) cos(m(/)) is a small perturbation 
with generally time-dependent amplitude R™(t) for the 
wave number m, which is integer in this case. Perturba- 
tion of the disk disturbs the magnetic potential as 

^{tj, t t) = ^-0(77, + *r (^7, t) cos(m0) + • • • , 

where 4*™ is the perturbation amplitude, which is associ- 
ated with the wave number to, and satisfies the equation 



^ Ml 7^) 
977! drj 




(21) 




z 



-1.5 -1 



Figure 9: Isolines of the base magnetic potential (r > 0) 
and the corresponding magnetic flux lines (r < 0) around a 
circular disk. 



Perturbation of the magnetic potential is related to that 
of the radius by the edge condition — resulting 



= -R\' 



dr 



Hi Uo 
Ro 1] 



(22) 



This perturbation vanishes with the distance from the 



(20) disk, i.e., 



0. Although Eq. (EU admits 



variable separation, such a solution is complicated by 
the edge singularity Ij22|l . Nevertheless, a compact an- 
alytic solution can be found with the following original 
approach. First, note that if ^' is a solution of the Laplace 
equation and e is a constant vector, (e • V)$ is a solution 
too. Second, if ^' satisfies a uniform boundary condi- 
tion Ijl6|l and e is directed along the boundary, (e • V)<i> 
satisfies that boundary condition too. Third, the op- 
erator (e • V) changes the radial dependence of 4" from 
~ (r — i?o)" to ~ (r — i?o)"~^, while the azimuthal depen- 
dence is changed from mode to to ?7i-|-l. Algebra becomes 
particularly simple when e is defined in the complex form 
as e — Ex + iCy = e"^(er + ie^). Then each application 
of the complex operator (e • V)<i> is accompanied by the 
multiplication with e"^. Thus, the solution for m = 1 
follows simply from the axisymmetric basic state as 



e-'^Ci (e • V) *o 



CiCq 



rj^ 



1/2 



Similarly, higher azimuthal modes can be found as 4"™ 



e ■ V) where is an axisymmetric solu- 
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tion satisfying Eq. H18|l . From the edge condition 

1 



we obtain 4'™ 



r/^™ ^ as ?7 



perturbation far away from the disk ^f" 
phes that along the disk '^q 



0. Moreover, vanishing of 
— > im- 



is a constant. The corresponding axisymmetric solution 
of Eq. Ijl8|l can be represented as 



k=l 

where Pn{x) and Qn[x) are the Legendre polynomials 
and functions of the second kind, respectively the 
expansion coefficients are c™ — {Ak — l)/™/(52fe-i(0), 
where 



T'ni 

— 



0F2i-2m(2m- 1)! 
(to - fc)!r(TO + fc + 1/2)' 



Then the solution for perturbation amplitude can be 
written as 



Dm-lD 



using the operator 



1 



m-2 ■ ■ 



(23) 



which is a spectral analog of (e • V) acting on the az- 
imuthal mode to. Calculation of Eq. lf23ll is algebrai cally 
complicated but can be done using mathematica 
which requires a considerable amount computer memory 
and, thus, is possible only for to < 5. Nevertheless, this 
suffices to deduce the general solution for arbitrary to. 



1 



m/2 



where the unknown constant Cm = BJ^ / Rq follows from 
Eq. 1I22I1 . It can easily be checked that the above solution 
indeed satisfies both Eq. II21II and the edge condition Ij22|l 
as well as the impermeability condition Ijl9|l . As for the 
semi-infinite sheet, the solution relative to the perturbed 
edge is obtained by the coordinate transformation 



*(r) =*(r' + e,.i?7') 



*o(r-') + ^^^^r + *r(0 

or 

*o(t-') + ^"(^')cos(to<?!)). 



where = '^Ti'^') - ^K*"')- I" the vicinity of 

the edge, this reduces to ^(7?,^) = C'o^y + 0[i^\ where 



Co = Co 



1- i(TO- l)i?7'/i?o 



. Note that there is no 



perturbation of the magnetic field with respect to the 
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Figure 10: Isolines of perturbation amplitudes of the mag- 
netic potential relative to the perturbed edge for the az- 
imuthal modes m = 2 (a) and m = 3 (b) plotted with the 
step 0.03 for C„Co = 1. 



edge for m = 1, because this mode corresponds to the 
offset of the disk as a whole. In this case, the field distri- 
bution moves together with the disk causing perturbation 
with respect to the original position of the disk, but not 
with respect to the disk itself. Perturbation amplitudes 
^™(r') are plotted in Fig. E3|for modes m — 2 and 3. 

The time-averaged force on the edge follows from Eq. 
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0, 



F = Fo + Fi^Fo[l-im^l) 



Ro 



where Fq — 7rCo/(8^o-Ro) = -So^o/(27r^o)- Similarly to 
the semi-infinite sheet, the normal stress balance at the 
edge of disk ifTTHl results in 



'dt 



r=Ro 



(24) 



where 7 is the surface tension of the disk. The hydrody- 
namic potential governed by Eq. ijHJ is found in cylindri- 
cal coordinates as $(r, (/), i) — $„i(t)r™ cos(m(/)), while 



the kinematic constraint Vr = d^/dr 



-Ro 



= dRX'/dt 



yields $™(t) l/{mR^-'^)dR'{'' /dt. The curvature per- 
turbation of the edge is 



rig 



rig 



(m^ — 1) cos(m(/)). 



Searching the edge perturbation as R^{t) — i?ie" 
where Am is in general a complex growth rate of the az- 
imuthal mode to, and substituting $ and Ki into Eq. 
Il24p . we eventually obtain 



to(to — 1) 



Bm 

TT 



(25) 



where tq — y/ pR^/j is the characteristic time of capil- 
lary oscillations; Bm = BqRq/{2p,qIqj) is the dimension- 
less magnetic Bond number characterizing the ratio of 
electromagnetic and surface tension forces. Without the 
magnetic field (Bm = 0), the growth rates for all modes 
are purely imaginary. Am — ±i(m — 1)\/toj correspond- 
ing to capillary oscillations of an inviscid disk. Increasing 
the magnetic field results in a decrease of the frequency 
of oscillations until the critical value of Bm is attained, 
at which an exponentially growing mode appears. Ac- 
cording to Eq. Ij25|l . the critical Bond number for mode 
TO, which is defined by the condition Am(Bm™) = 0, is 
Bm™ (m + l)7r, to = 2,3, . . .. Note that for to = 
and 1, we have Am — regardless of Bm because the first 
mode is not permitted by the incompressibility constraint 
for the layer of fixed thickness under consideration here. 
The mode to = 1, as already noted above, corresponds 
to the off'set of the disk as a whole, which has no effect 
relative to the disk itself as long as the external magnetic 
field is uniform. Thus, the first unstable mode is to = 2, 
for which the instability threshold is Bmc = Stt. Simi- 
larly to the straight edge case considered above, when 
Bm > Bmc the growth rate attains a maximum at the 
wave number TOmax defined by A^„^^^ = X'^^^_^^_i that 
yields TOmax = [|Bm/7r] , where the square brackets de- 
note the integer part. 




Figure 11: Top view of a flat Gallium drop in a transverse 
ac magnetic field of 13 kHz frequency at various induction 
amplitudes Bq: 12.9 mT (a), 23.5 mT (b), and 48.5 mT (c) . 



B. Comparison with experiment 

In the experiment, detailed description of which may 
be found in 01, a flat circular gallium drop of thick- 
ness do ~ 6 mm and radius -Ro ~ 30 mm was placed 
on a glass plate, which was slightly concave to center 
the drop, and put into a 6- winding solenoidal coil sup- 
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plied by ac current of / « 13 kHz frequency. At low 
currents in the coil, the drop was observed to be nearly 
circular, as seen in Fig. [TTT al. and remained such until 
the current exceeded some critical value, after which the 
drop became noticeably distorted, as seen in Fig. Illf bl. 
Further current increase resulted in the development of 
more corrugated drop shapes shown in Fig. Etc). Ac- 
cording to the experimental observations, the circular 
shape became unstable about the magnetic field induc- 
tion amplitude in the range 12.9 mT < Bq < 23.5 mT. 
Assuming the layer has approximately uniform curva- 
ture over the edge thickness with the radius ro ~ do/ 2 
that corresponds to an arclength ~ 7rdo/2, we find 
Bo,c = y/nBuicfJ.o'-fdo / Rq « 13.4 mT for the critical Bond 
number Bmc — 3n, where 7 = 0.718 N/m is the sur- 
face tension of gallium |i2(]j. This critical field strength is 
slightly higher than that in Fig. Ilir a) but considerably 
lower than that in Fig. HU h) . For the latter case we have 
Em « 30, which corresponds to the critical wave number 
rric — [Bm/n] — 1 = 8 defining the range of linearly unsta- 
ble modes 2 < m < rric- Note that the shape seen in Fig. 
[TTT b') has TO w 8, which corresponds to the critical mode 
for the given Bond number, although the fastest growing 
mode in this case is TOmax = 6. Given the simplicity of our 
theoretical model, these results may be thought to agree 
well with the experiment. There may be several reasons 
that preclude a better agreement with the experiment. 
First, the drop seen in Fig. Illf bl has a significant pertur- 
bation amplitude implying that its shape may be affected 
by nonlinear effects that are not accounted for by this lin- 
ear stability analysis. Second, our model of a thin per- 
fectly conducting sheet may be too simple for the given 
experiment with the relative drop thickness do/Ro « 0.2 
and the skin depth 6 — l/y/irfcrfio « 2.3 mm, where 
a = 3.7 X lO^il^^m^^ is the electrical conductivity of 
gaUium 

Note that our theory is developed for a disk of fixed 
thickness which excludes mode to = 0, while in the ex- 
periment the upper surface of the layer is free and this 
mode is permitted. Nevertheless, the theory is applica- 
ble also to this case because small-amplitude modes with 
TO > 1 are not coupled with the mode to = 0. The only 
difference is that the thickness of the layer may vary de- 
pending on the magnetic field. But once the equihbrium 
thickness is known, our theory can be used to predict 
whether the droplet will remain circular on the further 
increase of the magnetic field. 



IV. SUMMARY AND CONCLUSIONS 

We have analyzed the linear stability of the edge of 
a thin liquid metal layer subject to a transverse high- 
frequency ac magnetic field. The metal layer was consid- 
ered in the perfect conductor approximation supposing 
the ac frequency to be high so that the magnetic field is 
effectively expelled from the layer, while the thickness of 
the layer was assumed to be small relative to its lateral 



extension so that the layer could effectively be modeled 
as a thin perfectly conducting liquid sheet. First, we 
considered a general model of a semi-infinite sheet with 
a straight edge. This model, admitting an analytic solu- 
tion, allowed us to identify a pinch-type instability of the 
edge with the following simple mechanism. The magnetic 
field bending around the edge of a perfectly conducting 
layer creates a magnetic pressure on the edge trying to 
compress the layer laterally. In the basic state with a 
straight edge, the magnetic pressure, which is uniform 
along the edge, is balanced by a constant hydrostatic 
pressure in the layer. Perturbation of the edge in the 
form of a wave causes divergence of magnetic fiux lines 
at the wave crests and convergence in the troughs. This 
redistribution is because the magnetic fiux lines along the 
sheet are perpendicular to the current lines. But since the 
latter are aligned along the edge, the magnetic field has to 
be perpendicular to it. Consequently, the magnetic pres- 
sure is reduced at the crests and increased at the troughs, 
which drives the instability. Note that in this model of 
a thin sheet, the induction varies with the distance r 
from the edge as ~ l/v^ and, thus, it formally becomes 
singular at the edge. We circumvented this singularity 
by considering the sheet to have a small but neverthe- 
less finite thickness do. Then integration of the magnetic 
pressure, which scales as ^ l/c?o, over the thickness do 
resulted in a finite integral force on the edge indepen- 
dent of its actual thickness. This allowed us to obtain 
an analytical solution showing that the long-wave per- 
turbations are unstable when the wave number exceeds 
some critical value kc = Fo/{jIq), which is determined 
by the linear density of the electromagnetic force Fq act- 
ing on the edge, the surface tension 7, and the effective 
edge arclength ^o- The perturbations with wavelength 
shorter than critical are stabilized by the surface ten- 
sion, whereas the growth rate of long-wave perturbations 
reduces as ~ /c for fc ^ 0. Thus, there is the fastest grow- 
ing perturbation with the wave number /cmax = 2/3kc. 
When the layer is arranged vertically, long-wave pertur- 
bations are stabilized by the gravity, and the critical per- 
turbati on is c haracterized by the capillary wave number 
— \/ 9P/l- In this case, the critical linear density of 
electromagnetic force is Fq c = ^kJoj, which corresponds 
to the critical current amplitude /o,c = 4-y/7rfcc^o-^7/Mo 
when the magnetic field is generated by a straight wire 
at the distance L directly above the edge. Next, we solved 
analytically the linear stability problem for a thin circu- 
lar disk placed in a transverse uniform high-frequency ac 
magnetic field. It was found that the circular shape of 
the disk becomes unstable with respect to exponentially 
growing perturbation with the azimuthal wave number 
TO = 2 at the critical magnetic Bond number Bm^ = 37r. 
For Bm > Bm^ the wave number of the fastest growing 
perturbation is TOmax = [2Bm/(37r)] . These theoretical 
results were found to be in reasonably good agreement 
with available experimental data. 
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